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Abstract 

This  paper  describes  how  to  optimize  energy  storage  devices  (ESDs)  by  maximizing  their  net  present  value  (NPV).  This  requires  both 
technical  and  economic  information.  The  relevant  technical  information  is  specified  in  concise  form  by  the  energy-power  relation  (Ragone- 
plot)  of  the  ESD  and  its  lifetime.  The  economic  information  is  given  in  terms  of  operation  costs  (energy  costs),  investment  costs,  and  the 
economic  benefit  created  by  operating  the  ESD.  The  NPV  is  expressed  as  a  function  of  variables  such  as  the  size  of  the  ESD.  An  appropriate 
choice  of  these  variables  maximizes  the  NPV.  If  the  benefit  is  given  in  terms  of  the  energy  supplied  by  the  ESD  and  the  lifetime  is  operation 
independent,  the  optimization  reduces  to  a  minimization  of  the  total  lifecycle  costs.  In  this  case  a  knowledge  of  the  economic  benefit,  which  is 
often  the  quantity  that  is  most  difficult  to  model,  is  not  required.  In  more  general  cases,  e.g.  if  the  lifetime  depends  on  optimization  variables,  a 
quantification  of  the  benefit  is  necessary  in  order  to  maximize  the  NPV.  We  illustrate  the  approach  with  various  stationary  and  mobile 
applications,  and  for  batteries  and  capacitors.  ©  2002  Elsevier  Science  B.V.  All  rights  reserved. 
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1.  Introduction 

To  optimize  the  design  of  an  energy  storage  device  (ESD) 
means  different  things  to  different  people.  In  this  article,  it 
refers  to  a  maximization  of  its  net  present  value  (NPV)  [  1  ] . 
The  NPV  is  the  value  of  all  positive  cashflows  (revenues)  less 
the  value  of  all  negative  cashflows  (costs)  both  discounted  to 
a  single  date  to  make  them  commensurable.  Mathematically, 
the  NPV  is  a  weighted  sum  of  functions  containing  technical 
information  (efficiency,  power,  mass,  etc.),  where  the  weights 
reflect  economic  importance  in  monetary  terms  (costs  or 
revenues).  The  optimization  task  requires,  thus,  business 
knowledge,  and  engineering  know-how  in  an  interdisciplin¬ 
ary  manner.  Depending  on  the  particular  energy  storage 
design  problem  to  be  solved  it  can  be  anything  between 
straightforward  and  akward  to  translate  technical  duties  into 
monetary  terms.  A  straightforward  example  is  trading  of 
electrical  energy,  but  what  would  be,  e.g.  the  dollar  value  of 
having  a  generously  sized  starter  battery  in  a  car?  The 
concept  of  NPV  is,  however,  general  enough  to  allow  the 
description  of  a  wide  variety  of  energy  storage  problems. 

In  this  paper,  we  use  a  particular  simplifying  and  mean¬ 
ingful  description  of  ESDs  for  optimization:  energy-power 
relations,  the  so  called  Ragone  plots.  Ragone  plots  are 
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curves  which  relate  the  power  (density)  of  an  ESD  to  the 
available  energy  (density).  Due  to  internal  losses,  the  energy 
available  for  use  is  in  general  less  than  the  energy  stored  in 
the  ESD.  The  question  arises:  what  is  the  optimum  operation 
point  along  the  Ragone  curve  for  a  given  application?  Is  it 
the  maximum  possible  power  point  (“fast  discharge”),  or 
the  maximum  possible  energy  point  (“slow  discharge”),  or 
what  compromise  should  be  made  between  the  two?  This 
question  can  only  be  answered  within  a  combined  technical 
and  economical  description.  Ragone  plots  have  so  far  been 
mainly  used  for  a  rough  comparison  of  energy  storage 
technologies  across  orders  of  magnitude  in  either  power 
or  energy  capability.  However,  with  sufficient  care  in  the 
definition  and  sufficient  accuracy  in  the  measurement  of 
Ragone  plots,  they  may  serve  as  a  realistic  conceptual  tool 
for  the  actual  design  of  energy  storage  units.  On  the  other 
hand,  NPV  and  Ragone  plots  are  both  abstract  concepts.  In 
some  long  established  energy  storage  problems,  they  some¬ 
times  likely  produce  only  expected  or  even  trivial  answers. 
In  less  known  engineering  terrain,  however,  where  experi¬ 
ence  and  intuition  is  not  yet  mature  enough,  the  concepts 
discussed  in  the  present  paper  can  serve  indeed  as  a  sound 
starting  point  for  a  deeper  analysis. 

The  paper  is  organized  as  follows.  The  next  section 
provides  the  nomenclature  and  the  general  expression  for 
the  NPV  of  an  ESD.  Section  3  reviews  some  basics  of  Ragone 
plots.  Section  4  explains  how  to  solve  for  the  optimum  size 
of  the  ESD.  The  remaining  part  of  the  paper  is  devoted  to 
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illustrative  examples  going  from  easier  examples  to  more 
complex  ones.  Throughout  this  paper,  we  consider  two  types 
of  energy-power  relations,  which  describe  a  battery  and  a 
capacitor.  The  most  simple  application  (Section  5)  is  an  ESD 
for  a  stationary  application  with  constant  power  demand. 
After  providing  some  numerical  examples  in  Section  6,  the 
effect  of  a  simple  time-dependent  power  demand  is  discussed 
(Section  7).  Then,  in  Section  8,  we  investigate  a  mobile 
application,  where  the  power  is  partly  used  to  carry  the  ESD 
itself.  For  these  three  examples  the  benefit  is  given  in  terms  of 
the  energy  supplied  by  the  ESD,  and  a  constant  lifetime  is 
assumed.  As  a  consequence,  the  NPV  maximization  is 
equivalent  to  a  minimization  of  the  total  lifecycle  costs.  In 
Section  9,  we  generalize  our  approach  to  an  example  where 
the  lifetime  does  depend  on  the  operation  point.  This  case 
requires  information  on  the  economic  benefit  of  operating  the 
ESD  if  the  NPV  is  to  be  maximized.  As  an  example,  we 
compare  in  Section  9,  the  NPV  with  the  return  on  investment 
(ROI).  In  Section  10,  we  finally  discuss  a  different  direction 
of  generalization,  namely  a  dependence  of  the  benefit  on  the 
operation  point.  This  is  a  rather  general  case  that  will  be 
illustrated  for  an  ESD  for  uninteruptible  power  supply  (UPS). 


2.  The  net  present  value  of  an  energy  storage  device 

In  this  section,  we  derive  an  expression  for  the  NPV,  the 
sum  of  all  discounted  cashflows.  Because  the  purpose  is  to 
optimize  the  design  by  maximizing  the  NPV,  the  considera¬ 
tions  will  be  restricted  to  those  parts  that  depend  on  the 
design  variables  under  consideration.  The  absolute  NPV  is 
irrelevant  in  this  context  because  constant,  i.e.  design  inde¬ 
pendent  value  components  cannot  have  an  impact  on  design 
decisions  (except  the  decision  whether  to  have  an  ESD  at 
all).  Those  constant  components  vanish  on  differentiation. 
The  NPV  consists  of  the  present  value  of  all  benefits  Z?lc 
(revenues)  minus  the  present  value  of  all  costs  Clc  during 
the  lifetime  of  the  ESD: 

NPV  =  BLq  —  CLC-  (1) 

The  lifecycle  costs,  Clc,  are  usually  arranged  into  two  main 
groups,  initial  investment  costs,  Cjnv,  at  time  t  =  0,  and 
operation  costs,  Cop,  that  occur  during  the  system  lifetime  t: 

Clc  =  Cjnv  +  Cop.  (2) 


where  C?  is  the  size  independent  part  of  the  investment 
costs.  The  NPV  analysis  can  also  be  conducted  for  the  case 
of  a  non-linear  relationship  between  investment  costs  and 
size  and  also  for  the  dependence  of  investment  costs  on 
design  parameters  other  than  just  the  size.  Investment  costs 
should  also  include  the  present  value  of  the  costs  for 
dismantle  and  disposal  at  the  end  of  lifetime. 

Again,  we  make  a  restriction  that  can  be  easily  lifted  and 
assume  that  operation  costs,  Cop,„  incur  each  year  (0  <  n  <  t; 
n,  x  =  1,2, ...)  at  the  same  rate  during  the  product  lifetime 
t.  In  order  to  add  the  yearly  contributions  Cop„,  one  has 
to  calculate  their  present  value  at  the  time  of  investment. 
Denoting  by  r  the  cost  of  capital  in  percent  per  year,  one  has 


r  = 

'-^op 


E 


0  +  r)‘ 


(4) 


In  the  most  general  case,  both,  operation  costs  and  operation 
benefits  can  depend  on  a  variety  of  parameters.  For  the  sake 
of  clarity,  we  make  some  further  simplifying  restrictions.  We 
assume  that  the  operation  costs  can  be  separated  in  a  part 
associated  with  the  cost  of  the  energy  stored  in  the  ESD  and 
a  part,  Cj)p,  containing  the  rest  which  is  independent  of  the 
amount  of  energy  stored  (maintenance  costs,  etc.).  If  the  fuel 
costs  and  the  load  characteristics  are  the  same  for  each  year, 
and  if  aging  (i.e.  time  dependence)  of  the  efficiency  can  be 
neglected,  one  can  write  the  operation  costs  in  the  form: 

Cop  =  C°p  +  ced  [  ^  dr,  (5) 

where  Preq(t)  is  the  power  demand,  i.e.  the  power 
([Preq]  =  W)  delivered  by  the  ESD  to  the  load  at  time  f, 
ce  is  the  cost  of  energy  ([ce]  =  $/J),  and  i/  is  the  total 
(“round-trip”)  energy  efficiency,  i.e.  the  amount  of  energy 
that  is  delivered  from  the  ESD  divided  by  the  amount  of 
energy  used  to  charge  up  the  ESD.  This  means  that  the  total 
energy  efficiency  includes  the  charging  efficiency.  While  ;; 
depends  on  the  required  power  and  the  size  of  the  ESD,  the 
frequency  dependence  of  the  efficiency  is  neglected.  The 
case  where  the  time  variation  of  Preq(r)  is  so  fast  that 
frequency  dependence  becomes  important  will  be  discussed 
elsewhere.  Furthermore,  for  constant  r,  d  is 


5(i 


n=  1  (1  +  r) 


n—  1 


VT  -  1 
V-  1  ’ 


(6) 


The  unit  of  the  costs  will  be  [C]  =  $. 

We  restrict  our  analysis  to  ESDs  which  are  scalable  in  size. 
The  investment  costs  therefore  contain  a  term  dependent  on 
the  system  size.  We  furthermore  assume  this  dependency  to 
be  linear  such  that  the  ESD  consists  of  a  number  N  of  units 
with  costs  c,v  $  per  unit.  Examples  for  units  may  be  volume 
([N]  =  m3,  [cN]  =  $/m3),  mass  ([N]  =  kg,  [c^]  =  $/kg),  or 
number  of  cells  in  a  battery  or  capacitor  bank  ([(V]  =  1, 
[c/v]  =  $  per  cell),  etc.  The  investment  costs  are  then 

Cinv  =  C°v  +  CNN  (3) 


with  v  =  1/(1 -l-r).  Note  that  d  =  d( t)  is  a  function  of  the 
lifetime  and  may  thus  depend  on  the  operation  conditions.  In 
general,  this  fact  must  be  taken  into  account  in  an  optimiza¬ 
tion.  But  the  specific  dependence  is  rather  difficult  to  model, 
since  it  depends  in  general  on  the  full  load-profile  history. 
Only  in  the  case  of  a  constant  power  demand,  t  can  become 
a  simple  function  of  the  power  P.  Summarizing,  the  total 
lifecycle  costs  are 

Clc  =  C^  +  c^N  +  ced  f  — —  dr, 

J  1  year  ^7 


(7) 
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where  C\\c  contains  all  contributions,  which  are  not  con¬ 
tained  in  the  size  dependent  investment  costs  or  in  the  energy 
related  operation  costs.  For  simplicity,  we  will  omit  C\\c ,  it 
does  not  add  conceptual  difficulties. 

We  must  also  find  an  expression  for  the  total  monetary 
benefits  B\  q  due  to  providing  energy  from  the  ESD.  This 
expression  can  look  very  different  depending  on  circum¬ 
stances,  it  can  pose  modelling  challenges  being,  e.g.  a 
volatile  market  price  or  the  benefit  of  riding  through  a  power 
outage.  For  the  sake  of  clarity,  we  use  a  simple  description 
with  fixed  price  per  output  energy,  say  be  ([fic]  =  $/J).  Note 
that  the  pure  energy  margin  (i.e.  excluded  investment  costs) 
is  not  be  —  ce  but  be  —  ce/t],  since  there  is  always  an  energy 
loss  0 1  <  1).  With  assumptions  analogous  to  those  used  for 
the  derivation  of  the  operation  costs,  we  obtain 

BlC  =  M  f  Preq  d t.  (8) 

J 1  year 

The  remaining  task  of  this  paper  is  mainly  to  find  the 
maximum  of  Eq.  (1)  with  Eqs.  (7)  and  (8).  Before  we 
conclude  this  section,  we  mention  that  the  only  dependence 
of  /fix:  on  the  optimization  variables  in  our  restricted  sce¬ 
nario  enters  via  the  lifetime  through  d( t).  Flence,  if  the 
dependence  of  the  lifetime  on  the  operation  conditions  is 
negligible  (in  the  region  of  interest),  the  maximization  of  the 
NPV  reduces  to  the  minimization  of  Eq.  (7).  In  other  words, 
NPV  analysis  contains  the  lifecycle  cost  analysis  as  a  special 
case.  This  is  particularly  advantageous,  since  it  eliminates 
the  difficulty  of  modelling  BLq-  In  the  last  example,  we 
consider  the  general  case  of  an  operation  dependent  PL c, 
which  reflects  another  benefit  structure  of  the  ESD  for  a  user 
than  the  revenues  from  energy  supply. 


3.  Ragone  plots 

The  physical  information  on  an  ESD  comes  in  a  variety  of 
disguises.  It  is  connected  to  quantities  like  capacitance  or 
capacity,  internal  resistance,  efficiencies,  power  density, 
lifetime,  etc.  This  makes  a  comparison  between  energy 
storage  technologies  often  cumbersome.  Energy-power 
relations,  Ragone  plots,  have  a  two-fold  advantage  as  a 
concept  for  ESD  optimization:  they  are  rigorously  defined 
for  any  kind  of  ESD  [2]  and  they  readily  display  the  two 
parameters  with  cost  impact.  The  energy  efficiency  affects 
the  energy  costs,  while  the  investment  costs  for  a  specific 
application  depend  on  both  the  power  and  the  energy 
density.  The  information  on  efficiency,  energy,  and  power 
density  is  contained  in  the  Ragone  plot,  which  is  briefly 
reviewed  in  this  section. 

Energy  storage  technologies  are  characterized  by  a  typical 
power  scale  and  a  typical  energy  scale,  which  together 
determine  their  natural  field  of  applications.  For  example, 
batteries  have  high  energy  densities  and  are  used  for  long¬ 
time  applications  in  the  range  of  hours,  while  electric  double 
layer  supercapacitors  have  high  power  densities  and  are  used 


for  short-time  applications  down  to  fractions  of  a  second. 
It  is  convenient  to  compare  different  technologies  in  the 
energy-power  plane.  A  specific  ESD  does  not  correspond  to 
a  single  point  in  this  plane,  but  is  represented  by  a  curve, 
which  displays  the  energy  E  available  to  a  load  as  a  function  of 
the  power  P,  the  rate  with  which  the  energy  is  supplied  to  the 
load.  This  energy-power  relation  is  known  as  Ragone-plot 
[2-6].  The  exact  definition  of  an  energy-power  relation  may 
differ  among  different  publications  [7].  Usually,  these  rela¬ 
tions  are  plotted  on  a  log-log  scale  with  units  specific  energy 
([£]  =  J/kg)  versus  specific  power  ([P]  =  W/kg),  or  energy 
density  ([£]  =  J/m3)  versus  power  density  ([P]  =  W/m3). 

The  energy-power  relation,  E(P),  can  be  used  to  express 
the  energy  costs  in  Eq.  (7).  Here,  we  define  E  and  P  with 
respect  to  the  above  introduced  “unit”  of  N  (e.g.  specific 
energy,  energy  density,  energy  per  unit  cell,  etc.),  such  that 
the  power  P  delivered  by  a  single  unit  in  a  ESD  of  size  N  is 
PrCq//V.  If  Eq  is  the  initial  energy  (per  unit)  of  the  charged 
ESD,  E(P)/Eq  is  the  discharge  energy  efficiency.  Through¬ 
out  this  paper,  we  make  the  simplification  that  the  efficiency 
(in  contrast  to  the  lifetime)  does  not  depend  on  the  discharge 
depth.  By  introducing  a  charging  efficiency  i]c,  the  total 
energy  efficiency  can  finally  be  expressed  by  the  energy- 
power  relation: 


>1  = 


E(Preq/N) 

- F - '7c 


l(Preq/N). 


(9) 


Later,  we  will  assume  a  constant  charging  efficiency  >]c 
without  explicit  power  dependence,  except  for  a  single  exam¬ 
ple  with  charge-discharge  symmetry,  where  ;;c  =  E(P)/Eq, 
i.e.  charging  and  discharging  is  done  at  the  same  power. 

We  will  always  make  use  of  two  examples  of  ESD  for 
illustration  of  the  theory,  namely  the  (rechargeable)  battery 
and  the  (super)capacitor.  Both  are  characterized  by  their 
energy  capacity,  Eq ,  per  unit  and  by  their  maximum  power 
P max  per  unit.  The  internal  loss  will  be  described  by  an 
ohmic  resistance  P,  such  that  Pmax  =  Uq/AR,  where  Uq  is 
the  cell  voltage  of  the  charged-up  state.  The  energy  capa¬ 
cities  are  Eq  =  QqUq  for  the  battery  with  charge  capacity  Qq 
per  cell,  and  CI/q/2  for  the  capacitor  with  capacitance  C  per 
cell.  The  Ragone  plots  of  these  examples  are  derived  and 
discussed  in  Ref.  [2],  It  is  convenient  to  represent  them  in 
dimensionless  units,  p  =  P/Pmax  and  e{p)  =  E(pPmm)  /Eq, 
with  0  <  p,  e  <  1.  The  efficiency  as  a  function  of  p  is  given 
by  i]  =  rjce(p).  The  energy-power  relation  for  an  ideal 
battery  is  [2]  given  as 


e(p) 


1  P 

2(1 


(10) 


and  a  convenient  approximation  of  the  Ragone-plot  for  a 
capacitor-like  ESD  is  [7]  (we  do  not  consider  the  exact,  more 
complicated  expression  given  in  [2]) 

e(p)  =  l~P-  (ID 


The  two  functions  (10)  and  (11)  are  shown  in  Fig.  1.  The 
main  difference  between  them  is  that  the  available  energy 
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Fig.  1.  Normalized  Ragone  plots  for  a  battery  (solid  curve)  and  a  capacitor 
(dashed  curve).  The  energy  e  is  normalized  to  the  energy  capacity  (e  has 
then  the  meaning  of  a  discharge  efficiency),  and  the  power  p  to  the 
maximum  power.  The  dots  indicate  the  optimum  operation  point  for  K  —  1 
(cf.  Section  5). 


mentioned  problem  of  modelling  an  operation  dependent 
lifetime  t  as  a  function  of  the  whole  history  of  the  ESD,  we 
assume  that  it  can  be  modelled  by  a  reasonable  function  t (N) 
for  a  given  Preq.  Then,  Eq.  (13)  is  a  well-defined  function 
of  N  and  the  determination  of  the  maximum  is  straightfor¬ 
ward.  The  most  simple  case  refers  to  a  constant  lifetime  t 
independent  of  operation  power.  Optimization  of  the  NPV 
is  then  equivalent  to  a  minimization  of  the  lifecycle  costs. 
Combination  of  the  Eqs.  (12)  and  (13)  gives 


N 2 


cA(  ^j'(w«)d,  =  0 

Cn  J 1  year  T{Prcq/N) 


(14) 


where  the  prime  means  differentiation  with  respect  to  the 
argument  of  the  primed  function,  and  where  i]  is  given  by 
Eq.  (9).  One  observes  that  besides  the  efficiency  and  the 
load  function  Preq(?)  there  is  only  one  single  parameter 
that  characterizes  the  optimum.  In  dimensionless  units,  this 
parameter  has  the  meaning  of  the  ratio  of  energy  costs  to 
investment  costs.  For  more  general  cases,  additional  para¬ 
meters  may  enter.  In  the  following  sections,  specific  exam¬ 
ples  will  be  investigated. 


vanishes  at  maximum  power  for  the  capacitor,  while  the 
available  energy  remains  finite  for  battery-like  ESD;  more 
exactly,  it  is  one  half  of  the  full  energy.  Leakage  [2]  will  be 
neglected.  Leakage  leads  to  a  downturn  of  the  Ragone  curve 
at  very  low  power.  The  optimum  operation  point  of  a 
reasonable  ESD  is  usually  far  out  of  this  leakage  region. 
We  re-emphazise  that  our  examples  are  for  illustration  only; 
for  realistic  cases  the  optimization  should  be  carried  out  with 
measured  energy  efficiencies  or  measured  energy-power 
relations. 


4.  The  optimization  problem 


In  order  to  determine  the  optimum  design  of  an  ESD,  we 
express  the  NPV  as  a  function  of  the  design  variables. 
Consider,  e.g.  the  size  N  of  the  ESD.  Optimization  means 
then  replacement  of  P  by  PIcq  / N  in  the  NPV  and 

NPV(V)  =  max!  (12) 

where 


NPV  (V)  =  bed 


Pieq  df  —  Ced  [ 

J 1  year  *7 


df  —  cnN—  ClC. 


(13) 


The  first  term  on  the  right-hand  side  refers  to  the  benefit, 
which  is  related  to  the  energy  supplied  by  the  ESD.  This 
term  can  be  easily  generalized  and  can  be  an  arbitrary 
function  of  N  and,  if  necessary,  of  other  variables  appearing 
in  the  problem.  In  many  realistic  cases,  there  exists  a  local 
quadratic  maximum  (d(NPV)/d/V  =  0)  which  turns  out  to 
be  also  the  global  maximum.  In  order  to  avoid  the  above 


5.  ESD  for  constant  power  demand 


Consider  a  stationary  application  with  a  time  independent 
power  request  Preq  during  the  discharge  time.  The  required 
energy  per  year  is  aPreq  with  a  being  the  utilization  time 
per  year.  During  part  of  the  remaining  time  of  the  year, 
the  ESD  is  charged  with  a  charging  efficiency  r]c.  Using 
N  =  Prcq/pP  max,  the  NPV  can  be  written  in  the  form 


NPV  = 


Cv^req 


ej  p 


where  we  introduced  the  parameters 

CeCtdP  max 


K  = 
and 
P  = 


CnP  c 


Me 


(15) 


(16) 


(17) 


Note  that  the  units  of  Preq  and  Pmax  are  watt  and  watt  per  unit, 
respectively,  and  that  the  optimization  problem  is  indepen¬ 
dent  of  Preq.  The  meaning  of  K  is  the  ratio  between  total 
effective  energy  costs  per  unit  at  maximum  power  (more 
exactly:  the  lifetime  sum  of  energy  costs  for  the  specific 
utilization  factor)  and  the  investment  costs  per  unit.  The 
meaning  of  p  is  the  ratio  between  the  energy  benefit  (the 
energy  sales  price)  and  energy  costs  (taking  into  account  the 
charging  loss).  Since  d  is  constant,  maximum  NPV  corre¬ 
sponds  to  minimum  lifecycle  costs  which  can  be  written  as 


1 

Kp 


1 

e(p) 


min! 


(18) 
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The  only  parameter  occuring  in  this  optimization  problem 
is  K.  If  e(p)  is  a  decreasing  function  a  local  minimum 
of  Eq.  (18)  may  be  expected.  Minimization  of  Eq.  (18) 
leads  to 

Kp2e'  +  e2  =  0,  (19) 

where  the  prime  denotes  differentiation  with  respect  to  p. 
From  Eq.  (19)  one  finds  the  relation  between  K  and  the 
optimum  operation  point  p: 


p  =  i  (2  +  K-1  -  V4 K-1  +  K-2)  V4K-1  +  K~2,  (20) 


for  the  battery  (Eq.  (10)),  and 


P  = 


1 

T +7P 


(21) 


for  the  capacitor  (Eq.  (11)).  The  functions  p{K)  and  e(p(K)) 
are  plotted  in  Figs.  2  and  3,  respectively.  For  every  K  there  is 
an  optimum  operation  point  along  the  Ragone-plot.  The 
optimum  operation  point  for  K  =  1  is  indicated  in  Fig.  1  by 
solid  dots.  An  ESD  is  only  useful  if  the  maximum  NPV 
is  positive.  This  sets  a  lower  limit  to  the  energy  benefit 
(the  energy  price  obtained  by  discharging  the  ESD),  the 
minimum  energy  benefit  £>™ln.  The  relevant  relation  between 
/(  and  K  is  defined  by  the  equations  NPV  =  0,  and 
d(NPV) /d N  =  0.  The  solution  of  these  equations  must  in 
general  be  calculated  numerically.  Specific  examples  with 
the  assumption  that  Eq.  (15)  is  the  total  NPV  and  other  fixed 
costs  are  negligible  are  shown  in  Fig.  4.  One  observes,  that 
for  fixed  investment  costs  per  power,  the  minimum  energy 
benefit  for  batteries  is  lower  than  for  supercaps.  But 
note  that  usually,  supercapacitors  have  lower  cjv/Pmax- 
In  order  to  get  an  impression  on  the  typical  qualitative 
dependence  of  /?"""  on  the  various  parameters,  we  provide 


K 

Fig.  2.  Power  value  (p )  of  the  optimum  operation  point  as  a  function  of  K, 
for  constant  power  request  (battery:  solid  curve;  capacitor:  dashed  curve). 


K 

Fig.  3.  Energy  value  (e)  of  the  optimum  operation  point  as  a  function  of  K, 
for  constant  power  request  (battery:  solid  curve;  capacitor:  dashed  curve). 


the  result  for  the  capacitor,  which  can  be  analytically 
expressed: 


The  two  limit  cases  K  — >  0  and  K  — >  oc  are  as  one  expects. 
If  investment  costs  strongly  dominate,  the  present  value  of 
the  revenues  (per  unit)  have  to  be  at  least  equal  to  the 
investment  costs  (per  unit),  driPma,b"""  =  cn ■  On  the  other 
hand,  if  the  energy  costs  dominate,  the  energy  benefit  must  at 


Fig.  4.  Minimum  total  energy  price  as  a  function  of  the  investment  costs  per 
Watt  for  energy  costs  ce  —  0.1  $/kWh,  discount  rate  r  =  0.1,  and  utilizaton 
time  ot  —  2000  h  per  year.  Solid  curve:  lead-acid  battery  ( rjc  —  0.85, 
i  =  2.5  years).  Dashed  curve:  capacitor  ( rjc  =  0.95,  i  =  2  years). 
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least  equal  the  costs  of  the  stored  energy  h™n  =  ce/f7c 
(taking  into  account  the  charging  loss). 

Up  to  now,  the  considerations  were  restricted  to  a  given 
constant  charging  efficiency  i]c.  We  briefly  discuss  the 
important  case  of  charging-discharging  symmetry,  where 
;/c  =  e(p)  is  determined  by  optimization.  Using  this  in 
Eq.  (15),  optimization  leads  to 

2  Kp2e'  +  e3  =  0,  (23) 

where 

Ci-ocd 

K  =  Kr,  c  =  f—cN.  (24) 

*  max 

Eq.  (23)  should  be  compared  with  Eq.  (19)  and  leads  to 
qualitative  similar  results.  The  quantity  which  corresponds 
to  K  is  in  this  case  2K/e  or  2 Krjc/e.  In  other  words,  the 
operation  costs  have  two  times  the  weight  as  compared  to  the 
case  where  the  charging  efficiency  is  a  fixed  external  para¬ 
meter  that  is  not  influenced  by  the  design. 


6.  Numerical  examples 

Before  generalizing  the  previous  example,  we  provide  a 
few  typical  values  of  K,  the  optimum  operation  point,  and 
the  minimum  energy  benefit.  In  Table  1  various  battery  types 
are  listed  together  with  a  typical  supercapacitor.  Note  that 
there  are  many  different  types  of  lead-acid,  Ni-Cd,  and 
Ni-MeH  batteries  and  of  supercapacitors,  and  the  numbers 
presented  are  only  for  illustration.  We  assume  for  all  exam¬ 
ples  a  =  2000  h  utilization  time  per  year,  energy  costs 
ce  =  0.1  $/kWh,  and  an  interest  rate  r  =  10%.  We  consider 
Ni-MeH  and  Ni-Cd  batteries  which  are  similar  and  differ 
mainly  in  the  lifetime  and  the  efficiency.  It  turns  out  that  the 
lead-acid  battery  corresponds  to  the  highest  K  value  among 
the  batteries,  implying  that  it  has  to  be  operated  at  a  higher 
efficiency  than  the  Ni-Cd  or  the  Ni-MeH  batteries.  This  is 
mainly  due  to  the  low  investment  costs  of  the  lead-acid  cell. 
As  a  consequence,  it  has  the  lowest  required  energy  benefit 
among  the  batteries  in  the  table.  Of  course,  this  can  change 
for  mobile  (i.e.  transport)  applications  where  the  mass  plays 
an  important  role. 

In  real  life,  a  number  of  further  important  criteria  would 
play  decisive  roles  in  storage-device  selection:  self-discharge 
rates,  climatic  conditions,  safety  and  environmental  concerns, 


Table  1 

Comparison  of  optimization  results  for  different  ESD 


Lead-acid 

Ni-Cd 

Ni-MeH 

Supercapacitor 

Am*  (W/kg) 

100 

200 

200 

3000 

cN  ($/kg) 

5 

30 

30 

50 

t  (years) 

2.5 

4 

5 

2 

1c 

0.85 

0.8 

0.85 

0.95 

K 

11 

5.8 

6.5 

24 

P  =  pPmax 

45.2 

112 

108 

1022 

t>”m  ($/kWh) 

0.16 

0.19 

0.17 

0.13 

and  even  commercial  boundary  conditions  like  established 
supplier  relationships,  to  mention  a  few.  Although  this  level 
of  detail  is  beyond  the  scope  of  the  present  paper,  note  that 
most  of  these  criteria  can,  in  one  way  or  the  other,  be  taken 
into  account  in  the  NPV  analysis.  Self-discharge  (leakage) 
does  enter  the  Ragone-plot  [2],  or  it  can  be  taken  into 
account  in  an  effective  charging  efficiency.  Climatic  con¬ 
ditions  might  affect  the  lifetime.  Safety  and  environmental 
concerns  may  come  into  play  twice,  in  the  expression  for  the 
benefit  via  a  reduced  demand  and  in  the  expression  for  the 
total  lifecycle  costs  via  disposal  costs  at  the  end  of  life.  In 
any  case,  an  economic  comparison  of  different  ESD  types 
requires  the  complete  NPV,  including  also  constant  terms 
that  could  be  disregarded  for  an  optimization  of  the  single 
ESD. 

The  numbers  for  the  supercapacitor  should  not  be  com¬ 
pared  directly  with  the  batteries,  since  the  discharge  time  for 
supercaps  is  much  shorter,  and  thus,  the  two  different  types 
of  ESD  must  be  associated  with  completely  different  appli¬ 
cations.  The  assumption  of  identical  utilization  time 
a  =  2000  h  per  year  means  that  the  supercapacitor  having 
the  higher  power  capability  is  charged  and  discharged  with 
much  higher  frequency.  It  therefore  needs  a  lower  minimal 
energy  benefit.  Note  that  the  parameters  used  in  Fig.  4 
correspond  to  the  supercapacitor  and  the  lead-acid  battery 
of  Table  1. 


7.  ESD  with  time-dependent  power  demand 

If  the  power  demand  is  time-dependent,  Preq(t),  Eq.  (14) 
has  to  be  solved  numerically  for  most  realistic  cases.  We 
consider  a  specific  example  that  is  simple  enough  for  an 
analytical  treatment.  Assume  that  the  power  request  (per 
year)  is  equal  to  a  constant  P\  during  the  time  oq,  and  equal 
to  a  constant  P2  for  the  duration  <x2-  The  utilization  time  is 
now  a  =  oq  +  oq.  Without  restriction  we  take  P2  >  P\. 
Eq.  (14),  which  determines  the  optimum  size  N  of  the 
ESD  becomes 


K~l 


e'{yk ) 

e2(yk )  ’ 


(25) 


where  K  is  defined  in  Eq.  (16),  yic=Nic/N  with 
Nk  =  Pk/Pmax ,  and  cik  =  a*/a,  with  k  =  1  and  2.  It  is  clear 
that  N  >  N2  >  N\ .  The  optimum  system  size  N  as  a  function 
of  K  can  be  obtained  by  evaluating  the  right-hand  side  of 
Eq.  (25)  as  a  function  of  N  and  plotting  N  versus  K.  In  Fig.  5. 
We  plotted  the  system  size  N  for  the  battery  as  a  function  of 
K  for  various  values  of  eq  (=  1  —  a2),  and  for  N2  =  100, 
N\  =  10.  Cases  (a)  and  (e)  correspond  to  cases  with  constant 
demands,  P2  =  100  and  P\  =  10,  respectively.  There  is  a 
cross-over  between  (a)  and  (e)  for  N  >  100  at  high  K. 
However,  for  low  K  values  (high  investment  costs)  the 
optimum  system  size  is  limited  by  the  minimum  size  that 
can  supply  the  largest  power  request. 
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Fig.  5.  Optimum  system  size  as  a  function  of  K  for  an  ESD  that  has  to 
deliver  two  different  power  values,  P\  —  10Pmax  and  P2  —  100Pmax-  The 
different  curves  correspond  to  different  ratios  of  utilization  times 
ai(=  ai/(ai  +  0L2)).  (a)  ai  =  0,  (b)  ai  =  0.5,  (c)  oc\  —  0.9,  (d)  oq  =  0.999, 
(e)  oq  =  1.  For  small  K,  the  size  is  always  limited  by  the  maximum  power 
request  (N  >  100  except  for  (e)  where  N  >  10).  Besides  this,  there  is  a 
clear  cross-over  from  (a)  to  (e). 


K 

Fig.  6.  Optimum  operation  point  p  for  a  mobile  application  (solid  curve: 
battery;  dashed  curve:  capacitor).  The  ratio  of  power  needed  for 
transportation  of  the  ESD  unit  and  the  maximum  power  of  a  unit  is  assumed 
to  be  q  —  0.1.  The  thin  curves  correspond  to  the  case  q  —  0  for  comparison. 


with 


8.  ESD  for  transportation 


If  an  ESD  is  used  for  transportation  purposes,  one  has  to 
take  into  account  the  part  of  the  power  that  is  consumed  for 
carrying  the  ESD  itself.  This  implies  that  an  increase  of  the 
ESD  size  (weight)  leads  not  only  to  an  increase  in  energy 
storage  capability  but  also  to  an  increase  of  the  total  power 
demand.  For  example,  the  results  shown  in  Table  1  were 
established  for  stationary  applications;  it  is,  however,  well- 
known  that  lead-acid  batteries  are  an  inferior  technology  for 
electric  vehicles  due  to  their  heavy  weight.  We  denote  by 
P/eq  the  power  demand  of  the  vehicle  without  the  ESD’s 
weight,  and  write  the  total  power  demand  in  the  form: 


Preq  =  P°teq(l  +  KN),  (26) 

where  we  assume  a  linear  dependence  of  the  additional 
power  on  the  size  of  the  ESD.  The  quantity  zcP®  is  the 
additional  power  needed  for  transportation  of  one  unit  of  the 
ESD.  For  simplicity,  we  consider  again  a  constant  vehicle 
power  demand  P/eq  during  the  utilization  time  a.  The  rela¬ 
tion  N  =  Pleq/P  implies  that 


N  = 


p-*nq 


(27) 


In  the  dimensionless  units  (p  =  P/Pmm,  etc.),  the  optimiza¬ 
tion  condition  analogous  to  Eq.  (18)  becomes 


1 

p-q 


min! 


(28) 


q  = 


<eq 


(29) 


which  is  the  ratio  of  the  additional  power  needed  for  the 
transportation  of  the  ESD  unit  and  the  maximum  power  of 
the  unit.  Clearly,  it  must  hold  0  <  q  <  1  and  p  >  q,  because 
the  ESD  must  be  able  to  supply  the  power  for  carrying  its 
own  mass.  Solving  Eq.  (28)  leads  to 


e2(p) 

(q-p)pe'(p)  —  qe(p)  ' 


(30) 


While  this  function  can  probably  not  be  easily  inverted,  a 
plot  of  p  versus  K  is  shown  in  Fig.  6,  for  the  case  q  =  0.1. 
The  thin  curves  belong  to  q  =  0  and  are  identical  to  those 
in  Fig.  2.  As  one  expects,  the  solutions  for  mobile  and 
stationary  applications  differ  only  if  the  energy  costs  matter 
(K  — >  oo).  For  strongly  dominating  investment  costs 
(K  — >  0),  the  power  values  converge  towards  p  —>  1,  i.e. 
maximum  power  capability. 

The  optimum  power  p  for  dominating  energy  costs 
(K  — >  oo)  is  obtained  with  a  brief  calculation  from  Eq.  (30). 
Vanishing  of  the  denominator  requires  p  =  2^J~q  —  q  for 
the  battery,  and  p  =  v/cy  for  the  capacitor  model.  Using 
N  =  q/K(p  —  q)  from  Eq.  (27)  and  (29)  implies  that 

b  Jq  b  Ff 

N  = _ ^7 _ = _ rn.  (31) 

K  (1  -  y/q)  y/qO-y/ft  P max  ’ 

with  b  =  1/2  for  the  battery,  and  b  1  for  the  capacitor.  The 
limit  behavior  of  N  as  a  function  of  q  for  dominating  energy 
costs  can  easily  be  understood;  q  =  0  is  identical  to  the  case 
of  a  stationary  application  and  shows  the  same  behaviour 
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(N  — »  oo  for  K  — »  oo ).  On  the  other  hand,  q  —  I  means  that 
the  power  remaining  for  transportation  of  the  bare  vehicle 
(without  ESD)  goes  to  zero,  leading  again  to  IV  — >  oo. 

We  mention  that  there  is  a  variety  of  further  interesting 
ESD  optimization  problems  in  transportation.  For  example, 
instead  of  modelling  the  payback  per  energy  unit,  the  sales 
price  that  can  be  obtained  by  a  customer  could  be  a  function 
of  the  maximum  vehicle  range  and  of  available  acceleration. 
Then,  also  be  would  be  a  function  of  p  and  can  be  included  in 
the  optimization  in  a  straightforward  way.  This  would, 
however,  go  beyond  the  scope  of  this  paper. 


9.  ESD  with  an  operation  dependent  lifetime 

Up  to  this  point,  the  lifetime  used  in  our  considerations 
was  a  fixed  parameter  determined  by  other  components  in 
the  system  than  the  ESD.  In  this  section,  we  consider  the 
case  where  the  lifetime  refers  to  the  time-to-breakdown  of 
the  ESD  itself.  Now,  z  and,  consequently,  d  and  K  depend  on 
the  operation  power  p. 

Considering  a  constant  power  Preq  during  the  utilization 
time  a,  a  maximization  of  the  NPV  (15)  yields 

K^j  e2  =  0  (32) 

with  K'  =  dK /dp.  We  assume  a  lifetime  model  that  leads  to 
K  =  K{)(  1  —  p)‘ ,  where  the  ESD  breaks  immediately  (i.e. 
the  lifetime  vanishes)  if  the  power  reaches  Pmax;  Kq  and  y  are 
model  parameters.  We  assume  a  constant  discharge  depth 
for  all  applications,  such  that  the  only  reason  of  a  lifetime 
decrease  is  due  to  an  increased  power.  Reasonable  ESD  have 
y-values  with  0  <  y  1.  The  NPV  as  a  function  of  the 
operation  point  is  plotted  in  Fig.  7,  for  y  =  0.5  and  in  units 


Kp1  e'  +  (  1  +  p2  (  p 


Fig.  7.  Net  present  value  as  a  function  of  p  for  the  case  where  the  lifetime 
depends  on  p.  Parameter  values:  Kq  —  1,  p  —  10;  solid  curve:  y  —  0.5; 
dotted  curve:  y  —  0  (constant  lifetime);  dashed-dotted  curve:  ROI 
(dimensionless  units)  for  y  =  0.5. 


of  CjvPreq/^max-  It  is  compared  to  the  case  with  constant 
lifetime  (dotted  curve).  The  NPV  is  lower  and  the  maximum 
is  shifted  to  lower  values  of  the  operation  power,  as  one 
expects  because  higher  power  now  implies  shorter  lifetime 
in  addition  to  lower  efficiency. 

Let  us  now  compare  maximization  of  the  NPV  with  a 
maximization  of  the  return  on  investment  (ROI),  defined  as 
the  ratio  of  the  NPV  to  the  present  value  of  all  costs, 
NPV/Clc-  The  ROI  is  shown  by  the  dashed-dotted  curve 
in  Fig.  7.  A  straightforward  maximization  yields: 

Kp2e'+  (l  +plpje2  =  Q,  (33) 

from  which  a  relation  between  the  optimum  operation  point 
p  and  Kq  is  obtained.  The  optimum  operation  point  as  a 
function  of  Kq  is  shown  in  Fig.  8,  for  NPV  optimization 
(solid  curve)  and  ROI  optimization  (dashed-dotted  line).  We 
assumed  for  illustration  p  =  10.  Note  that  in  general  there 
exists  a  lower  bound  for  the  NPV  or  ROI  in  the  sense  that 
values  below  this  bound  are  not  acceptable,  since  the 
investment  costs  become  too  high.  This  bound  is  not  neces¬ 
sarily  at  zero,  since  there  might  be  other  contributions  to  the 
NPV,  which  are  independent  of  the  energy.  Nevertheless,  we 
indicate  by  a  dot  in  Fig.  8  the  point  where  the  NPV  and  the 
ROI  vanish.  From  the  definitions  of  NPV  and  ROI  it  follows 
directly  that  this  point  is  at  the  same  time  the  intersection 
point  of  theses  curves.  The  result  for  constant  lifetime,  being 
equivalent  to  the  result  for  minimization  of  lifecycle  costs,  is 
also  indicated  for  illustration  (dotted  line).  The  operation 
point  optimizing  the  NPV  is  shifted  to  lower  values,  as  one 
expects.  The  ROI  approach  leads  to  results  that  differ  only  at 
low  Kq  significantly  from  the  approach  of  minimization  of 


Fig.  8.  Optimum  operation  point  p  as  a  function  of  Kq  with  y  —  0.5, 
p  —  10,  and  for  the  case  where  the  lifetime  depends  on  p  (solid  curve). 
Dashed-dotted  curve:  operation  point  from  ROI  maximation;  the  black  dot 
indicates  the  zero  of  the  NPV  and  ROI.  Dotted  curve:  y  —  0  (constant 
lifetime). 
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lifecycle  costs,  i.e.  at  large  power.  This  is  understood  by 
recognizing  that  in  our  case  this  approach  can  be  interpreted 
as  a  minimization  of  average  costs  per  effective  lifetime, 
C\  c/d  =  minimum.  At  small  p  this  approximates  well  the 
constant  lifetime  case.  At  p  close  to  1,  the  decrease  of  d  leads 
to  an  increase  of  C^c/d.  Consequently,  the  optimum  opera¬ 
tion  point  is  significantly  lower  than  for  the  constant  lifetime 
case.  The  largest  possible  value  of  p  for  maximum  ROI  is 
given  by  1/(1  +  y).  There  is  generally  a  large  difference 
between  the  result  of  the  NPV  optimization  and  the  ROI 
optimization.  We  conclude  this  section  with  the  remark,  that 
for  p  <2  the  example  shown  in  the  figure  (y  =  1/2)  can 
have  a  global  non-quadratic  maximum  of  the  NPV  which  is 
at  the  boundary  p  =  1.  Since  this  example  is  mainly  for 
illustration  we  will  not  go  further  into  details. 

10.  ESD  with  an  operation  dependent  benefit 


Fig.  9.  Optimum  operation  point  p  as  a  function  of  r);  battery  (solid  curve) 
and  a  capacitor  (dashed  curve).  The  dots  indicate  the  vanishing  of  the  NPV, 
which  becomes  nageative  at  low  <5. 


Until  now  the  benefit  Z?lc  has  been  associated  with  the 
monetary  value  of  the  energy  supplied  by  the  ESD.  Benefits 
can,  however,  be  given  in  terms  of  other  values  that  character¬ 
ize  the  utility  derived  by  using  the  ESD.  In  this  section,  we 
model  the  benefit  by  the  price  a  customer  is  willing  to  pay  for 
the  ESD.  This  price  depends  on  the  operational  performance. 

For  illustration,  consider  an  ESD  used  as  an  uninterrup¬ 
tible  power  supply  (UPS).  The  UPS  is  obviously  more  useful 
the  longer  its  discharge  time  (for  constant  Pre q).  As  short 
power  interruptions  are  much  more  frequent  than  long  ones, 
the  price  of  the  ESD  increases  sublinearly  with  the  discharge 
time,  t.  We  model  the  dependence  of  the  benefit  on  the 
discharge  time  by  Blc(0  =  BoPleqy/t  with  a  constant  Bq. 
The  energy  costs  are  negligible  for  an  UPS  such  that  the 
restriction  to  K  =  0  is  appropriate.  Then,  the  NPV  becomes 


NPV  = 


C.V  Prcq 


p  p 


(34) 


where  <5  =  BoPmax  yj (£o/Pmax)  / cv  is  the  ratio  between 
the  customer  benefit  per  unit  for  the  characteristic  dis¬ 
charge  time  £o/Pmax,  and  the  investment  costs  per  unit. 


It  holds  8  =  Bq \J (QqUq/4R) / cn  for  the  battery  and 
<5  =  Bq\J (C / AR)  Uq/cn  for  the  capacitor  (cf.  Section  3). 
The  NPV  is  easily  optimized  leading  to  a  relation  between 
optimum  operation  point  and  the  parameter  <5: 


8  = 


'I'feJp 

e  —  pe' 


(35) 


The  relation  (35)  is  plotted  in  Fig.  9,  for  both  battery  and 
capacitor.  For  the  capacitor,  the  simple  expression 
p  =  4/(4  +  S2)  is  obtained.  Obviously,  for  the  same  8  the 
capacitor  is  optimal  at  lower  p  than  the  battery  due  to  the 
lower  energy  efficiency  which  leads  to  a  shorter  discharge 
time.  In  the  figure,  the  points  where  the  NPV  vanish  are 
indicated  by  dots.  The  values  of  (p.  8)  corresponding  to  zero 
NPV  are  (0.889, 1.3)  and  (0.5,2)  for  the  battery  and  the 


capacitor,  respectively.  The  absolute  value  of  the  NPV  will, 
in  general,  contain  constant  contributions  which  are  not 
included  in  Eq.  (34). 

We  conclude  this  section  by  providing  a  numerical  example. 
Consider  an  ESD  with  E0  =  30  Wh/kg,  Pmax  =  100  W/kg, 
Cjv  =  5  $/kg.  A  typical  value  for  Bq  is  250$/kWv/h.  One 
obtains  then  8  =  2.1  with  a  corresponding  optimum  opera¬ 
tion  point  p  =  0.45.  In  other  words,  the  optimum  mass  of  the 
battery  is  22  kg. 


11.  Conclusion 

We  have  shown  that  the  energy-power  relation  of  an  ESD 
(Ragone-plot)  can  be  used  for  the  calculation  of  the  optimum 
operation  point,  or  the  optimum  size  of  an  ESD  for  a  given 
application.  A  number  of  illustrative  examples  have  been 
discussed.  In  the  most  simple  case,  the  only  economical 
parameter  appearing  is  the  ratio  between  the  energy  costs 
and  investment  costs.  This  parameter  determines  the  balance 
between  maximum  efficiency  and  maximum  power.  The 
optimization  depends  on  the  specific  application.  We  illu¬ 
strated  this  for  applications  with  dynamic  load  and  for  mobile 
applications.  For  fixed  lifetime  and  for  a  benefit  given  in 
terms  of  the  energy  supplied  by  the  ESD,  the  optimization 
task  is  equivalent  to  a  life  cycle  cost  minimization.  On  the 
other  hand,  a  lifetime  which  depends  significantly  on  the 
operation  point  requires  knowledge  on  the  benefits  in  order  to 
optimize  the  NPV.  It  has  been  shown,  that  the  optimization  of 
the  ROI,  which  does  not  require  knowledge  on  the  revenues, 
may  lead  to  results  that  differ  significantly  from  the  NPV 
optimum.  For  optimization  of  a  commercial  product,  the 
NPV  optimization  should  be  preferred  from  a  rational  point 
of  view.  Finally,  we  discussed  a  simple  case  of  an  UPS  where 
the  benefit  depends  on  the  operation  point. 

It  should  be  clear  that  the  present  approach  is  rather  general, 
and  our  examples  can  be  generalized  in  a  straightforward  way. 
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For  example,  the  NPV  of  an  ESD  for  electric  vehicles  can 
now  easily  be  written  down,  with  a  benefit  depending  on  the 
range  (discharge  time)  and  the  accelleration  (related  to  the 
power),  and  including  operation  dependent  lifetime  as  well 
as  energy  costs.  One  can  also  consider  several  variables 
leading  to  a  multi-dimensional  optimization  problem,  and 
one  can  include  constraints.  This  is  from  a  conceptional 
point  of  view  straightforward,  and  leads  to  a  realistic 
optimum  design.  In  many  cases,  our  approach  can  easily 
be  extended  to  other  types  of  power  supplies.  For  example, 
fuel  cells  are  treated  similar  to  batteries  with  taking  for  >/  the 
total  system  efficiency  as  a  function  of  power  density. 
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